HARDY-LIEB-THIRRING INEQUALITIES FOR FRACTIONAL 
SCHRODINGER OPERATORS 



RUPERT L. FRANK, ELLIOTT H. LIEB, AND ROBERT SEIRINGER 

Abstract. We show that the Lieb-Thirring inequalities on moments of negative 
eigenvalues of Schrodinger-like operators remain true, with possibly different con- 
stants, when the critical Hardy-weight C|a;| -2 is subtracted from the Laplace oper- 
ator. We do so by first establishing a Sobolev inequality for such operators. Similar 
results are true for fractional powers of the Laplacian and the Hardy-weight and, 
in particular, for relativistic Schrodinger operators. We also allow for the inclusion 
of magnetic vector potentials. As an application, we extend, for the first time, the 
proof of stability of relativistic matter with magnetic fields all the way up to the 
critical value of the nuclear charge Za — 2/tt. 



Contents 



1. 



r^tochrcMorJ 



2. Main Results! 



2.1. Hardv-Lieb-Thirring inequalities! 

3^_Sj)bj3k2y^Ha£dj_^ 

3.2. Proof of Theorem Ipj 

3.3. Proof of Corollary [ 2. 5J 
4. Grou nd S^ate Representation! 



5^ i ^roof j of i thc - ljaKl^Ljx^ Inequalities! 

5.2. Heat kernel estimate ) 

5.3. Proof of Theorem I2~l1 



6, . Extension to Ma gnetic Fields! 
(^^^^A^r^to£d^m^qu;mt^ 

6.3. Extension of Theorem 2. II to magnetic fields! 



Appendix A. A Constant in the Sobolev Inequality (12.61) 


ADDendix B. Localization Formula in the Magnetic Case 


ADDendix C. Closabilitv of the Quadratic Form bn A 


ADDendix D. Nash's Argument 




ADDendix E. A Trace Estimate 


References! 





Date: October 26, 2006. 
© 2006 by the authors, 
purposes. 



10 
11 
12 
14 
15 
15 
16 
IS 
IS 
19 
21 
21 
21 
24 
25 
26 
27 
2S 



This paper may be reproduced, in its entirety, for non-commercial 



2 



RUPERT L. FRANK, ELLIOTT H. LIEB, AND ROBERT SEIRINGER 



1. Introduction 

We shall generalize several well known inequalities about the negative spectrum of 
Schrodinger-like operators on IR d . As an application of our results we shall give a proof 
of the 'stability of relativistic matter' — one which goes further than previous proofs 
by permitting the inclusion of magnetic fields for values of the nuclear charge all the 
way up to Za = 2 /it, which is the critical value in the absence of a field. 

There are three main inequalities to which we refer. The first is Hardy's inequality, 
whose classical form for d > 3 is the following. (In this introduction we shall not be 
precise about the space of functions in question, but will be precise later on.) 

f (\Vu{x)\ 2 - ( d ~^ 2 \u(x)A dx>0. (1.1) 
The second is Sobolev's inequality for d > 3, 

/ \Vu(x)\ 2 dx>S 2 , d / \u(x)\ 2d ^ d - 2 Ux\ = S 2 , d \\u\\j d/{d _ 2) . (1.2) 

The third is the Lieb-Thirring (LT) inequality |LThj for the Schrodinger operator 
H = —A — V(x). If its negative eigenvalues are denoted by — Ai < — A2 < • • • , and if 
7 > 0, then 

E A J^W V{x)l +d/2 dx = L 1 AV + \Ct- ( L3 ) 

j JR d 

This holds if and only if 7 > | when d — 1, 7 > when d = 2 and 7 > when d > 3. 
(Here and in the sequel t_ := max{0, — t} and t + := max{0,t} denote the negative 
and positive parts of t.) 

By duality, (jl.2j) is equivalent to the fact that the Schrodinger operator H = — A— V, 
with has no negative eigenvalues if ||V + ||d/2 < ^.d- On the other hand, ()1.3|) gives an 
upper bound to the number of negative eigenvalues in terms of ||V+||<2/2 when 7 = 
and it estimates the magnitude of these eigenvalues when 7 > 0. 

All three inequalities can be generalized by the inclusion of a magnetic vector po- 
tential A (related to the magnetic field B by B = cuxlA). That is, V is replaced 
by V — iA(x), and A by (V — iA(x)) 2 . The sharp constants in p.ljl . p. 2)1 remain 
unchanged while the sharp constants in p. 3)1 that are independent of A might, in prin- 
ciple, be different from the ones for A = 0. However, the best constants known so far do 
not depend on A. The inclusion of A is easily done in p.ljl . p. 2)1 by using the diamag- 
netic inequality, but the inclusion in p. 3)1 is more delicate; one uses the Feynman-Kac 
path integral formula to show that for each i,i/6l <! and t, r > 0, the A-field reduces 
the magnitude of the heat kernel e t< ^~ tA ^ (x, y) relative to e tA (x, y), and hence reduces 
the resolvent kernel |[— (V — iA) 2 + r] _1 (x,y)| relative to [—A + r] _1 (x,?/). 
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From |BVaj one can deduce that (jl.lj) and (|1.2jl can be combined as follows: For 

e > 0, 

J (W(x)| 2 - ^rHx)\ 2 + K*)| 2 ) dx > S' 2A£ \\u\\l d/(d _ 2+e) . (1.4) 

with S' 2 d e — >■ as e — > 0. Note the extra term |w(x)| 2 on the left side to account for 
the fact that the left and right sides behave differently under scaling; examples show 
that it really is necessary to have e > here. 

In |EkFr j . a parallel extension of (jl.3|) for the negative eigenvalues — Xj of the 
Schrodinger operator H = —A — (d — 2) 2 /(4|x| 2 ) — V is proved. For 7 > and 
d > 3, 

E A J ^ ^ I V(x)l +d/2 dx = L' l4 \\V + \\]Xm > ( L5 ) 

with L' d > L 1>d . Note that there is no need for an e in (|1.5|) : the fact that e ^ 
in ()1.4|) is reflected here in the fact that 7 > is needed. As before, a magnetic 
vector potential can easily be included in ()1.4|) . but it does not seem easy to include 
a magnetic field in ()1.5|) by the methods in |EkFr j . 

Our goal here is to extend these results in several ways. One extension is to include 
a magnetic field in (|1.5j) . Another is to consider fractional powers of the (magnetic) 
Laplacian, i.e., to the case in which |V — iA\ 2 is replaced by |V — iA\ 2s with < s < 
min{l, d/2} (which means that we can now include one- and two-dimensions). This is 
a significant generalization because the operator (— A) s is not a differential operator 
and it is not 'local'. Really different techniques will be needed. In particular, we shall 
use the heat kernel to prove the analog of (jl.5)l , in the manner of |Llj . A bound on this 
kernel, in turn, will be derived from a Sobolev-like inequality (the analogue of (|l-4jl ) 
by using an analogue of Nash's inequality, as explained in |LLoj . The appropriate 
inequalities are naturally formulated in a weighted space with measure |x| _/3 (ix for 
j3 > 0. Therefore, a pointwise bound on the heat kernel for a weighted 'Hardy' 
operator exp{— t\x\ a ((— A) s — C Std \x\~ 2s + l)|x| Q } for appropriate a > will be needed 
and will not be straightforward to obtain. 

In the dimension most relevant for physics, d = 3, the earlier case s — 1 may 
be called the non-relativistic case, while the new result for s — 1/2 may be called 
the relativistic case. Indeed, the resulting LT inequality, together with some of the 
methodology in |LYj . yields a new proof of the stability of relativistic matter, which 
will be sketched in Section l2~2l The main point, however, is that this new proof allows 
for an arbitrary magnetic vector potential A. Since the constant in the relativistic 
(s = 1/2) Hardy inequality that replaces (d — 2) 2 /4 is 2/n (which is the same as 
the critical value of Za in the field-free relativistic case), we conclude that we can 
simultaneously have an arbitrary A-field and the critical value of Za, the nuclear 
charge times the fine-structure constant. Up to now it was not possible to have both 
an arbitrary field and Za = 2/tt. Therefore, the proof of the analogue of ()1.5|) for 
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s = 1/2 with an A-field opens a slightly improved perspective on the interaction of 
matter and radiation. 

Acknowledgements. We thank Heinz Siedentop for suggesting that we study inequal- 
ities of this type, and we thank him, Ari Laptev and Jan Philip Solovej for helpful 
discussions. We also thank Renming Song for valuable comments on a previous version 
of this manuscript. This work was partially supported by the Swedish Foundation for 
International Cooperation in Research and Higher Education (STINT) (R.F.), by U.S. 
National Science Foundation grants PHY 01 39984 (E.L.) and PHY 03 53181 (R.S.), 
and by an A. P. Sloan Fellowship (R.S.). 

2. Main Results 

2.1. Hardy-Lieb-Thirring inequalities. We recall the Hardy-type inequality 

f \x\- 2s \u(x)\ 2 dx < C'l [ |£i 2s K0| 2 ^, ueC™(R d ), (2.1) 
valid for < 2s < d. Here 

:= (2ny d/2 [ u(x)e-^ x dx, £ E R d , 

denotes the Fourier transform of u. The sharp constant in 1)2. 1J) . 

_ T%d + 2sm 
M r 2 ((d-2s)/4) ' [ ] 

has been found independently by Herbst |Hej and Yafaev [Yaj. Moreover, in the latter 
paper it is shown that this constant is not achieved in the class of functions for which 
both sides are finite. In the case < s < min{l,cf/2}, this fact can also be deduced 
from our ground state representation in Proposition 14. 1[ which therefore represents an 
independent proof of (|2.1jl . See Remark 14.21 below. 

We denote D = — iV. Consider a magnetic vector potential A e L 2 QC (M. d , M. d ) and 
the self-adjoint operator (D — A) 2 in L 2 (W d ). For < s < 1 we define the operator 
\D — A\ 2s := ((D — A) 2 ) s by the spectral theorem. One form of the diamagnetic 
inequality states that if < s < 1 and u 6 dom\D — A\ s then \u\ G H s (M. d ) and 

||(-A) s/2 |m||| 2 < \\\D- A\ s u\\ 2 . (2.3) 

Here and in the sequel, || • || = || • |h denotes the L 2 -norm. We refer to Remark |6.2I for 
more details concerning ()2.3|) . Combining this inequality with the Hardy inequality 
(|2.1|) we find that the quadratic form 

h sA [u] := \\\D- A\ s u\^ -C s4 \\\x\- S u\\ 2 (2.4) 

is non-negative on dom|D — A\ s if < s < min{l,(i/2}. We use the same notation 
for its closure and denote by 

H sA = \D - A\ 2s - C s4 \x\- 2s 
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the corresponding self-adjoint operator in L 2 (R d ). 
Our main result is 

Theorem 2.1 (Hardy-Lieb-Thirring inequalities). Let 7 > and < s < 

min{l, d/2}. Then there exists a constant L jAs > such that for all V and A 

tr (\D - A\ 2s - C s4 \x\~ 2s -V)l< L lAs [ V{x)\ +d,2s dx. (2.5) 

Remark 2.2. For d > 3, Theorem 12. II holds also for s — 1. In the non-magnetic case 
A = 0, this has been proved in |EkFr j . In the present paper, we present an independent 
proof of this result, which allows for the inclusion of an A-field. For simplicity, we 
restrict our attention to s < 1 in the following, and comment on the (simpler) case 
s — 1 in Subsection 16.41 

Theorem 12.11 will be proved in Section |3J The main ingredient in its proof is a 
Sobolev-type inequality which might be of independent interest and which we shall 
present in the remainder of this subsection. 

For the case A = we shall drop the index A from the notation, i.e., 

h s [u}= [ |£| 2s K0| 2 ^-C M / \x\~ 2s \u(x)\ 2 dx. 

We denote the closure of this form by the same letter. In particular, its domain 
domh s is the closure of C£°(~R d ) with respect to the norm (h s [u] + HwH 2 ) 1 / 2 . Note that 
H s (M. d ) C domh s with strict inclusion. In particular, there exist functions u G dom/i s 
for which both sides of (j2.1)l are infinite. 

Hardy's inequality (|2.1j) implies that h s is non-negative. The following theorem 
shows that for functions of compact support it even satisfies a Sobolev-type inequality; 
i.e., h s [u] can be bounded from below by an L 9 -norm of u. 

Theorem 2.3 (Local Sobolev-Hardy inequality). Let < s < min{l,c?/2} and 

1 < q < 2* = 2d/(d — 2s). Then there exists a constant C q ^, s > such that for any 
domain Q C M. d with finite measure \Q\ one has 

\\u\\ 2 q < C qAs \Q\ 2 ^h s [u], u e C^(Q). (2.6) 

Remark 2.4. Note that the exponent q is strictly smaller than the critical Sobolev 
exponent 2* = 2d/(d— 2s). By considering functions which diverge like \x\~( d ~ 2s ^ 2 at 
x = 0, it is easy to see that the inequality (j2.6j) cannot hold if q is replaced by 2*. 

We note that the analogue of Theorem 12 .31 in the local case s — 1 is proved in |B Vaj . 
Theorem 12 . 31 will be proved in Sectional where we also deduce the following corollary 
from it. 

Corollary 2.5 (Global Sobolev-Hardy inequality). LetO < s < min{l,c?/2} and 

2 < q < 2* = 2dj id — 2s) . Then there exists a constant C d s > such that 

Ml < C' qA s(hs[u] + INI 2 ), u G L7 °°(R d ). (2.7) 
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Note that (J2.7)) may be written in the scale-invariant form 

NI^c^am^^ihiKI-^), (2.8) 

where C^ d s can be expressed explicitly in terms of C' q d s . This inequality follows from 
applying ()2.7|) to functions of the form U\(x) = u(Xx) and then optimizing over the 
choice of A. 

Although the sharp constants in Theorems 12.11 and 12.31 as well as in Corollary 12.51 
are unknown, explicit upper bounds involving a certain variational expression can be 
deduced from our proof. In Appendix El we do evaluate explicit bounds on the con- 
stants in Theorem 12 .31 in the special case d — 3, s — 1/2, which is the most interesting 
case from a physical point of view, as will be explained in the next subsection. 

Remark 2.6. Corollary 12. 51 is one of the main ingredients in our proof of Theorem l2.il 
On the other hand, Corollary 12.51 is an easy consequence of Theorem 12. 1| except for 
the value of the constants. In fact, ()2.5|) implies that 

h s [u]> I V(x)\u(x)\ 2 dx- (l iAs [ V(x)l +d/2s dx 

for all u e domh s and all V E L 1+d / 2s (M. d ) . Eq. (j2.8|) follows by optimizing the right 
side over all V. 

2.2. Stability of relativistic matter. We shall now explain how the inequalities in 
Theorem 12.11 can be used to prove stability of relativistic matter in the presence of 
an external magnetic field. The proof works up to and including the critical value 
Za = 2/tt, which is a new result and solves a problem that has been open for a long 
time. We refer to |L2[ IL3j for a review of this topic. 

We consider N electrons of mass m > with q spin states (q = 2 for real elec- 
trons) and K fixed nuclei with (distinct) coordinates Ri, . . . ,Rr £ and charges 
Zi,...,Z K > 0. A pseudo-relativistic description of the corresponding quantum- 
mechanical system is given by the Hamiltonian 

H n ,k = X) ( V (Pi ~ V® A ( x j)) 2 + m2 ~ m J + aV NtK (x 1 , ...,x N ). 
j=i V / 

The Pauli exclusion principle dictates that Hn,k ac ts on functions in the anti-symmetric 
iV-fold tensor product of L 2 (M 3 ; C 9 ). Here we use units where h = c = 1, a>0is the 
fine structure constant, and 

N K 

v NjK {xi, . . . ,x N ) := ^2 \ x i ~ ^r 1 ~ y. ) z k\ x j - Rki' 1 

l<i<j<N j=l k=l 

+ ^3 Z^ZilRk — Ri\ 1 . 

l<k<l<K 

Stability of matter means that Hn,k is bounded from below by a constant times 
(N + K), independently of the positions Rk of the nuclei. 
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By combining the methods in |LYj and our Theorem 12 .![ one can prove the following 

Theorem 2.7 (Stability of relativistic matter). There is an a c > such that for 
all N, K , qat < a c and a max{Zi, . . . , Zk} < 2/7r one has 

H N , K > -mN. 

The constant a c can be chosen independently of m, A and Ri, . . . , Rk- 

The constant a c in Theorem 12 . 71 depends on the optimal constant in the Hardy- LT 
inequality (j2.5|) for d = 3 and s = 1/2. A bound on this constant, in turn, can be 
obtained from our proof in terms of the constant in the Sobolev-Hardy inequalities 
(|2.6|) . We do derive a bound on the relevant constants in Appendix [21 but these 
bounds are probably far from optimal. In particular, the available constants do not 
yield realistic values of a c so far. 

After completing this work, we discovered a different proof of the special case of the 
Hardy- LT inequality needed in the proof of Theorem 12 .7\ namely the case where the 
potential V is constant inside the unit ball, and infinite outside. In this special case 
a substantially improved constant can be obtained, and this permits the conclusion 
that Theorem 12 . 71 holds for the physical value of a, which equals a ~ 1/137. We refer 
to |FLSj for details. 

We briefly outline the proof of Theorem 12.71 An examination of the proof in |LYj 
shows that there are two places that do not permit the inclusion of a magnetic vector 
potential A. These are Theorem 9 (Localization of kinetic energy - general form) and 
Theorem 11 (Lower bound to the short-range energy in a ball). Our Lemma lB.ll in 
Appendix[B]is precisely the extension of Theorem 9 to the magnetic case. This lemma 
implies that Theorem 10 in |LY| holds also in the magnetic case, without change except 
for replacing \D\ by \D — A\. 

Our Theorem 12.11 can be used instead of Theorem 11 in |LYj . In fact, the left side 
of (3.20) in |LY| is bounded from below by qilxIlL times the sum of the negative 
eigenvalues of \D — A\ — — Ci?~ 1 6'^(x), where Or denotes the characteristic 

function of a ball of radius R. By Theorem I2.1| this latter sum is bounded from 
below by a constant times C^R^ 1 . The resulting bound is of the same form as the 
right side of (3.20), except for the constant. It is this constant that determines the 
maximally allowed value of the fine structure constant, a c . The rest of the proof 
remains unchanged. Note that, in particular, the Daubechies inequality |Dauj remains 
true also in the presence of a magnetic field. 

2.3. Outline of the paper. Before giving the proofs of our main results, we pause 
to outline the structure of this paper. 

• In the next Section 01 we give the proof of the Sobolev-Hardy inequalities in 
Theorem 12.31 and Corollary 12.51 

• In the following Section|3]we prove what is customarily called the "ground state 
representation" in Proposition 14. 1[ except that here the "ground state" fails to 
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be an L 2 function. Such a representation for fractional differential operators 
does not seem to have appeared in the literature before. 

• In Section 03 we give the proof of our main Theorem 12.11 about Hardy-LT 
inequalities. We first consider the non-magnetic case A = 0. One of the key 
ingredients is the ground state representation obtained in Section HJ which 
allows us to prove a certain contraction property of the heat kernel in some 
weighted L 1 -spaces. Nash's argument |LLo[ Sect. 8.15-18] then allows us to 
translate the Sobolev-Hardy inequalities in Corollary 12. 51 into pointwise bounds 
on the heat kernel in an appropriate weighted L p -space. These bounds lead to 
the Hardy-LT inequalities via the trace formula in Proposition 15.31 in the spirit 

of |EH- 

• Finally, in Section |U] we derive diamagnetic inequalities which will allow us to 
extend the proof of Theorem 12. II to the magnetic case. 

3. Sobolev-Hardy Inequalities 

Our goal in this section is to prove Theorem 12.31 and Corollary 12.51 We start with 
a short outline of the structure of the proof. 

Our proof is based on the fact that we can control the singularity of H s ip near the 
origin if we know the singularity of ip at that point (cf. Lemma Theorem 12 . 31 fol- 
lows by observing that the L 9 -norm of a symmetric decreasing function can be bounded 
above by integrating the function against \x\ d ^ l / q ^ 1 \ see Lemma 13.41 Moreover, it is 
enough to restrict one's attention to symmetric decreasing functions. Corollary 12.51 
follows from Theorem 12.31 by an IMS-type localization argument, see Lemma [3.51 

We present some auxiliary results in the following Subsection 13.11 The next two 
Subsections 13.21 and 13.31 contain the proofs of Theorem 12.31 and Corollary 12.51 respec- 
tively. 

3.1. Auxiliary material. We start with the following integral representation of the 
operator (— A) s . 



Lemma 3.1. Let d > 1 and < s < 1. Then for all u G H 

1 2s 



|£H«(0| a de = fl.,d/ / W ^_^} 2 dxdy, (3.1) 



where 



n . o2,-i - rf / 2 r((d + 2s)/2) 
a s4 .= 2 vr | r( _ s) | ■ (3-2) 

Lemma 13.11 is well known; we sketch the proof for the sake of completeness. 

Proof. For fixed y we change coordinates z = x — y and apply Plancherel. Recalling 
that (u{- + z)Y(i) = e^ z u(£) we obtain 

\u{x) -u(y)\ 2 



\x-y 



d+2s 
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The integral in brackets is of the form c S! d|£| 2s , with 

/*oo p 

c s4 := / / \e iru - e -l\ 2 der- 28 - x dr 

Jo JS^ 1 

POD 

= 2 / (IS^ 1 ! - (2n) d/2 r- (d -^ 2 J {d _ 2)/2 (r)) r" 2 ^ 1 dr . 



Here, Jid-2)/2 is the Bessel function of the first kind of order (d — 2)/2 |AbStj . Recall 
that (S^ 1 ! = 2n d / 2 /T(d/2). The formula (JH2J) for c B)d = follows now from 

r- (J M/2 (r) - 2-(^)/ 2 r(d/2)-V^ 2 )/ 2 ) rfr = 2 -^^| 

for d/2 <Rez<(d + 4)/2, see jYa[ (2.20)]. ■ 

Let us recall that \x\~ a is a tempered distribution for < a < d with Fourier 
transform 

b a (I " r) A (0 = b d _ a \£\- d+a , b a := 2«/ 2 r(«/2) (3.3) 

(see, e.g., |LLo| Thm. 5.9], where another convention for the Fourier transform is 
used, however). We assume now that s < d/2. Then (—A) s \x\~ a is an L^-function 
for < a < d — 2s and 

((-A) s -C s4 \x\- 2s ) \x\' a = $ s4 (a)\x\- a - 2s , (3.4) 

where C s> d is defined in ()2.2|) and 



$ a Aa) := : C s . 





V2sbd-a 


b d - 


h 


2 2s 


(Y{{a 


V r((d 



(3.5) 



Later on we will need the following information about the a-dependence of 

Lemma 3.2. The function <& s d is negative and strictly increasing in (0, (d — 2s)/2) 
with $ s ,d((d - 2s)/2) = 0. 

Proof. First one checks that 

lim $ M (o;) = -C s>d < 0, $ Sjd ((d - 2s)/2) = 0. (3.6) 

a— >0 

Now we abbreviate (3 := a/2, r := d/2 and write 

/(/?) := T((3)/T(r - (3), g((3) := /(/3 + 

so that $ Si d(a:) = 2 2s g(/3) — C s , d - In view of ()3.6|) it suffices to verify that g{(3) is 
strictly increasing with respect to (3 € (0, (r — s)/2). One finds that 
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with ip = T'/r the Digamma function. Hence 

where, in view of ()3.7j) . h{t) := ip'(t) — i[)'(r — t). Since ip' is strictly decreasing (see 
[ATStl (6.4.1)]), one has h(t) > for t G (0,r/2). This proves that g\f3) > for all 
(3 G (0, (r — 2s)/2). In the case /3 G ((r — 2s)/2, (r — s)/2) one uses in addition the 
symmetry fo(t) = — /i(r — t). ■ 



3.2. Proof of Theorem 12.31 Our proof of Theorem 12.31 is close in spirit to [BL 
where a remainder term in the Sobolev inequality on bounded domains was found. 
We first exhibit functions ip G domh s in the form domain which do not lie in the 
operator domain but for which the singularity of the distribution (indeed, function) 
H a ip at x = can be calculated explicitly. 

Lemma 3.3. Let < x < 1 be a smooth function on M. + of compact support, with 
x(r) = 1 for r < 1. Define 

:=x(M/A)|a;|- Q (3.8) 

for < a < (d — 2s) /2 and A > 0. Taen ^ G domh s for < s < 1 and, for every 
e > 0, there exists a \ £ = X £ (a, d, s, x) such that for any A > A e; 

(((-A) s - C s4 \x\- 2s ) ipx) (x) < -(|$ M (a)| - e) \x\- a ~ 2s for all x G B, (3.9) 

in the sense of distributions. Here, $ Sit j(a) is given in h3. .5j) . and £> denotes the unit 
ball in R d . 

Proof. It is not difficult to show that ip\ G fP(!R d ), which implies that ip\ G dom/i s . 
(Consult the proof of Proposition 14. II for details.) Let < (p G C^°(B). According to 
(|3.4j) one has 

(V>a, ((-A) s - C M |x|" 2s ) y>) = $ M (a) (|x|- Q - 2s , <p) - (^ A , (-A) >' 

where ip\(x) := (1 — x(M/A))M~ Q - It follows from Lemma 13.11 (with the aid of 
polarization) that 

(&,(-A)V) = -2^// *^^>-p(A) jf,^*, 

with 

p(A) = su P 2a M |xr +2s / | My l 2s dy= sup 2a s4 \x\° +2s [ nsr%*- 

(3.10) 

Note that p(X) is finite for A > 1, and monotone decreasing to as A — > oo. Hence, 
for a given e > we can choose X e such that p(X e ) = e. Since is negative by 

Lemma I3~2l we have established (|3.9j) . ■ 
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Lemma 3.4. Let 1 < q < oo and u G L q (M. d ) a symmetric decreasing function. Then 



\u 



<q- l \B\- l/q ' u(x)\x\' d/q ' dx (3.11) 



where \B\ is the volume of the unit ball B in IR d , and 1/q + 1/q' = 1. 

Proof. First note that ()3.11j) is true (with equality) if u is the characteristic function of 
a centered ball. For general u we use the layer cake representation |LLo| Thm. 1.13], 
u(x) = f^°Xt(x)dt, where Xt is the characteristic function of a centered ball of a 
certain i-dependent radius. Then, by Minkowski's inequality |LLo[ Thm. 2.4], 



\ U \\q 



poo POO P 

< \\ X t\\ q dt = q- 1 \B\- l/q ' / X t(x)\x\- d/q ' dxdt 
Jo Jo J 



JO 

l \B\- l,q ' f u(x)\x\- d/q ' dx, 



= q 

proving ()3.11|) . ■ 
Now we give the 

Proof of Theorem \2.'A We remark first that we may assume Q to be a ball and u to 
be a spherically symmetric decreasing function. Indeed, passing to the symmetric 
decreasing rearrangement of u leaves the left side of ()2.(ij) invariant while it decreases 
the right side. The kinetic energy term on the right side is decreased by virtue of 
Riesz's rearrangement inequality (compare with |LLo[ Thm. 3.7, Lemma 7.17]) and 
J \u\ 2 \x\~ 2s dx increases |LLo[ Thm. 3.4]. Moreover, by scaling we may assume that 
Q = B, the unit ball. 

Since B is bounded, Holder's inequality implies that it suffices to prove ()2.6|) for 
d/(d — 2s) < q < 2d/(d — 2s). For such q let a := d/q' — 2s and note that < 
a < (d — 2s)/ 2. It follows from Lemmas 13.41 and 13.31 that for symmetric decreasing 
functions u on B 

\\u\\ q < g^lSr 179 ' / u(x)\x\- d/q ' dx < 2q- 1 \B\- 1/q '\<l? s4 (a)\- 1 \(u,H s ^)\. 
Jb 

Here ip = is chosen as in Lemma 13.31 with e = |$ Sj rf(a)|/2. An application of 
Schwarz's inequality, \(u, H s ip)\ 2 < h s [u}h s [ip], concludes the proof of (|2.6j) . ■ 

In Appendix |XJ we shall give an upper bound on the constant appearing in the 
Sobolev inequality in the special case of d = 3 and s = 1/2, which is the case of 
interest in the application in Section l2~2l 



3.3. Proof of Corollary 12. 5L We will deduce Corollary 12.51 from Theorem 12.31 bv a 
localization argument. For comparison, we recall first the IMS formula in the local case 
s = 1. If xo, ■ ■ ■ j Xn are Lipschitz continuous functions on M d satisfying Y^=o Xj = 1> 
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then 

„ l2 {d-2) 2 \u\ 2 \ . 

4 |x| 2 / 

= E / (|V(^)| 2 - (rf - 2)21 ^ 12 ) dx-f ±\v Xj \M 2 dx. 

j=Q JM d V 4 \ X \ J JR d j=0 

The following analogous formula for the non-local case, suggested by Michael Loss, is 
given in |LYj . The proof is an immediate consequence of Lemma 13.11 For a general- 
ization to the magnetic case, see Lemma [B. II below. 

Lemma 3.5. Let < s < min{l,cf/2} and let xo,---,Xn be Lipschitz continuous 
functions on M. d satisfying YTj=o x] = 1 ■ Then 

n 

h s [u] = J2h s [xM - (u,Lu), u e C™(R d ), (3.13) 

3=0 

where L is the bounded operator with integral kernel 

n 

L(x,y) := a s4 \x - y\- d ~ 2s J2(Xj(x) ~ Xj{y)) 2 - 

3=0 

Let us recall the following (non-critical) Sobolev embedding theorem, which is easy 
to prove. (Cf., e.g., the proof of |LLo| Thm. 8.5].) If s < d/2 and 2 < q < 2* = 
2d/(d - 2s) then H s (R d ) C L q {R d ) and 

\\u\\ 2 q < S qAs (||(-A) S / 2 M || 2 + || M || 2 ) , u G H s (R d ). (3.14) 

In combination with the localization Lemma [3.51 this allows us to give the 

Proof of Corollaru \2.5l Let Xo> Xi De smooth functions on R d with xl + Xi = 1 such 
that xo(x) = if |ar| > 1 and xi(x) = if |x| < 1/2. Let 2 < q < 2d/{d - 2s). Then, 
by Theorem 12.31 

||Xo«|i; < C qAs \B\ 2 (^)h s [ Xo u}, 

and by (j3~T4]l 

llxi«|i;< S qAs (U-Ay/ 2 (xiu)\\ 2 + \\xiu\\ 2 ) 

< S qAs (h s [ X iu] + (2 2s C s4 + i)||xm|| 2 ) . 

Hence Corollary 12 . 5l follows from Lemma E31 noting that L is a bounded operator. ■ 



4. Ground State Representation 

Eq. ()3.4|) and Lemma 13.21 suggest that the function \x\~( d ~ 2s ^ 2 is a 'generalized 
ground state' of the operator H s . Our next goal is to establish a ground state repre- 
sentation. Let us recall the analogous formula in the 'local' case s = 1. If d > 3 and 
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v(x) = \x\( d 2 ^ 2 u(x) then 

The corresponding formula in the non-local case < s < 1 is more complicated but 
close in spirit. It was derived some years ago by Michael Loss (unpublished notes) for 
the relativistic case s = 1/2 and d = 3. 

Proposition 4.1 (Ground State Representation). Let < s < min{l, d/2}. If 

u e Cg°(R d \ {0}) and v(x) = \x\^ d - 2s ^ 2 u(x) , then 

h\u]=a J f \< x )~<y)\ 2 _± d JL_ U 2) 

41 * k *h* \x-y\ d+2s \ X \(d-2s)/2 \y\(d-2s)/2 ' ^ 

with a s ^ given in 

Proof. Let < a < (d - 2s)/ 2. We shall prove that if u G C7^°(M d \ {0}) and 
v a (x) := |x| q m(x) then 



f \Z\ 2s \u{0\ 2 d£-(C Sid + $ s>d {a)) [ \x\- 2s \u(x)\ 2 dx 

JR d JR d 

\v a (x) - v a (y)\ 2 dx dy 



(4.3) 

\x - y\ d + 2s \x\ a \y\ a 

The proposition follows by letting a — *■ (d — 2s)/2. Indeed, the constant in front of 
the second integral on the left side then converges to C S}d , according to Lemma 13.21 
By splitting the integral into four regions according to the support of u, it is easy to 
see that the right side is continuous in a and converges to the right side of (j4.2|) . 

For the proof of (|4.3|) we can suppose that the support of u is in the unit ball. We 
shall first prove the equality for mollified versions of |x| _c \ namely functions w n (x) = 
\x\~ a x(x/n), where x £ C£°(WL d ) with x( x ) = 1 f° r \ x \ — 1- 

Let us first show that u n e H s (M. d ). It is clearly in L 2 (M. d ), hence it suffices to es- 
tablish that (— A) s / 2 uj n £ L 2 (R. d ). According to |LLo| Thm. 5.9] the Fourier transform 
of u n is given by the convolution of x an d Since x is assumed to be smooth, x 

decays faster than any power of |£|. It is then easy to see that u n decays like 
and hence € L 2 (R d ). 

By polarization in Lemma f3. 11 we get for any / and g in H s (M. d ), 



\e s M)m = a s , d 1 1 (f(x) - f(y))(g(x) - g{y)) -^%^- s . (4.4) 

<R d \ x — y\ 



We apply this formula to g(x) = u n (x) and f(x) = \u{x)\ 2 / uj n {x) = \u(x)\ 2 \x\ a . In 
this case, the right side of (|4.4jl is given by 



a s,d / / \u(x) - u(y)\' 



u(x) u{y) 



uj n (x) uj n (y) 



uj n (x)uj n (y) ) | ^^ +2s • (4.5) 
x y\ 
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Note that u{x)/uj n {x) = u(x)\x\ a = v a (x) is independent of n, and is a Cq° function 
since the origin is not in the support of u by assumption. By dominated convergence, 
(|4.5jl converges to 

a s> d ff (Hx) - u(y)\ 2 - \v a (x) - v a (y)\ 2 \x\~ a \y\~ a ) ^^ +2s (4.6) 
J JR d xR d \ x ~~ y\ 

as n — ► oo. The left side of (|4.4J) can be written as (compare with ()3.3j) ) 

(2vr) d / 2 ^ // iei 2s 7^n d x(^(e-n)irr- d ^^'. (4.7) 

Since / decays faster than polynomially, | ■ \ 2s f G L p (R d ) for any 1 < p < oo. Hence 
its convolution with the approximate 5-function (2ir) d / 2 n d x(n ■ ) converges to | ■ \ 2s f 
strongly in any L p , for 1 < p < oo |LLo| Thm. 2.16]. Therefore, ()4.7|) converges to 

b^a f |f|2 ^ d( = ^A_« r Hx)l 2 lxl -2s dx _ (4 g) 

Here we used (13. 3|) again. The equality of ()4.6|) and ()4.8|) proves f|4.3|) . ■ 

Remark 4.2. From the representation ()4.2|1 we immediately recover Hardy's inequal- 
ities ()2.1)1 in the case s < min{l,<i/2}. Moreover, we see that equality can not be 
attained (for a non-zero function). From ()4.2|) . it is also easy to see that the constant 
C s ^d is sharp. For this, consider a sequence of functions u n , supported in B, approx- 
imating | | — C rf — 2s )/ 2 close to the origin in a suitable sense. The right side of ()4.2|) 
remains finite in the limit n — ► oo, whereas J \u n (x)\ 2 \x\~ 2s dx diverges. 

5. Proof of the Hardy-Lieb-Thirring Inequalities 

This section contains the proof of our main result in Theorem 12. II We consider here 
only the non-magnetic case A = 0, the extension to non-zero A will be straightforward 
given the necessary diamagnetic inequalities which we derive in the next Section |U1 
We explain the necessary modifications in the proof of Theorem 12 .11 in Subsection 16.31 

The ground state representation (j4.2j) suggests that it is more natural to regard 
h s [u] as a function of v given by v(x) = \x\( d ~ 2s ^ 2 u(x). In terms of this function v, 
the Sobolev-Hardy inequality in ()2.7|) can be formulated in the weighted space with 
measure {x^^dx, where (3 = q(d — 2s) /2. Namely, 

IMli^Rd,!*!-^) ^ C 'q,dA V ^ B f 3V ^>L 2 (R d ,\x\-Pdx) , (5.1) 

where Bp is the operator on L 2 (R d , \x\~ ^dx) defined by the quadratic form 

{v,B0o)v mx \- Pdx) = h s [\x\^ d - 2s ^ 2 v} + \\\x\- {d - 2s)2 v\\ 2 L2{Rd4x) . (5.2) 

We suppress the dependence on s in Bp for simplicity. Note that the right side of (|5.2jl 
is independent of (3. The dependence of Bp on (3 comes from the measure \x\~^dx of the 
underlying L 2 space, which is determined by the value of q in the Sobolev inequality 
(|5-H) as (3 = q{d — 2s)/2. We emphasize again that the choice of the weight \x\^^ 2s ^ 2 
on the right side of ()5.2|) is determined by the ground state representation ()4.2|) . 
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The proof of Theorem 12.11 proceeds in the following steps. 

• From the ground state representation ([4.2)1 . we will deduce that Bp satis- 
fies the Beurling-Deny criteria, which implies that e~ tBl3 is a contraction on 
L x (M. d , \x\~Pdx) (Subsection El) . 

• Together with the Sobolev-Hardy inequality (J5.1)) this yields a bound on the 
kernel of e~ tBl3 via Nash's method (Subsection 15. 2|) . 

• This bound on the heat kernel can then be translated into a LT bound in the 
spirit of |Llj (Subsection 15. 3|) . 

5.1. Contraction property of Bp. Let Sjp := L 2 (R d , |x| _/3 <ir). We assume that 
d — 2s < (3 < d, which corresponds to 2 < q < 2* in f[5.1|) . The quadratic form 



considered in the Hilbert space Sjp, is non-negative and closable on C^°(R d \ {0}), and 
hence generates a self-adjoint operator Bp in ftp. 

We shall deduce some positivity properties of the operator exp(— tBp). By Propo- 
sition 14.11 the quadratic form bp satisfies 

(1) if v ,w G dombp are real-valued then bp[v + iw] = bp[v) + bp[w], 

(2) if v G dombp is real- valued then \v\ G domfr^ and bp[\v \) < bp[v], 

(3) if v G dombp is non-negative then min(i>, 1) G dombp and bp[min(v, 1)] < bp[v]. 

By a theorem of Beurling-Deny (see [Dj Section 1.3] or |ReSi2| Section XIII. 12]) this 
implies that exp(— tBp) is positivity-preserving and a contraction in L 1 (M d , \x\~^dx). 
That is, it maps non-negative functions into non-negative functions, and it decreases 
L 1 -norms. 

5.2. Heat kernel estimate. From the contraction property derived above, we will 
deduce a pointwise bound on the heat kernel, i.e., on the kernel of the integral operator 
exp(— tBp). We emphasize that this kernel is defined by 



(exjp(-tBp)v)(x) = / exp(-tBp)(x,y)v(y)-—rg . 

jR d \y\ 

We shall use the Sobolev inequality (|5.1|) for this purpose. 

Proposition 5.1. Let d — 2s < (3 < d. Then exp(—tBp) is an integral operator on 
ftp and its kernel satisfies 



Proof. Let 9 := (q — 2)/(q — 1) G (0, 1). Then Holder's inequality and Corollary 12.51 
yield for any v G C^{R d \ {0}) 



bp[v] ■=h s [\x\- {d - 2s)/2 v] + \\\x 



(d-2a)/2.,||2 




< \\V 



1-9 

L3(M d ,|x|-' 3 (ix) 



e 

L 1 (M d ,|x|-' 3 (ia;) 




,/(l-0)/2 



bp[v]^\\v\\ 



e 



L 2 (M d ,|a;| _ ' 3 dx) 
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Equivalently, if p is as in the proposition, then 

This is a Nash-type inequality in R d with measure \x\~@dx. By Nash's argument (see 
|LLo| Theorem 8.16] or jDj Section 2.4]) this implies that exp(— tBp) is an integral 
operator with kernel satisfying ()5.3|) . with the constant Kp^, s given in the proposition. 
For the sake of completeness we sketch the proof of this claim in Appendix O below. 



Remark 5.2. In the above argument we used the Nash-type inequality (J5.4j) which 
we had deduced from the Sobolev-type inequality (|2.7|) . The heat kernel bound (|5.3|) 
would actually follow directly from the latter inequality by [HI Thm. 2.4.2]. However, 
we preferred the simplicity of the above argument, yielding in addition an explicit 
constant. 

5.3. Proof of Theorem 12. It Step 1. As a first step, we seek an upper bound on 
the number of eigenvalues below — r of the operator H s — V, which we denote by 
N(—r,H s — V). By the variational principle we may assume that V > 0. Then the 
Birman-Schwinger principle (see, e.g., |ReSi2j ) implies that for any increasing non- 
negative function F on (0, oo) 

iV(-l, H s - V) < F(l)" 1 tr F (V 1/2 (H S + iy l V 1/2 ) . 

Let IA : L>2(W i ) — > Sjp be the unitary operator which maps u i— > |x|^ 2 u. Then 

V 1/2 (H S + /)- V 1/2 = WW X J 2 B^ X W X J 2 U , (5.5) 

where Wp is the multiplication operator on Sjp which multiplies by the function 
Wp(x) := \x\P +2s - d V(x). Therefore 

trF (W 2 {H S + J)-VV 2 ) = tr % F (w^B^W 1 / 2 ) . (5.6) 

We need the following trace estimate. 

Proposition 5.3. Let f be a non-negative convex function on [0, oo), growing poly- 
nomially at infinity and vanishing near the origin, and let 



DC 



F(X) := / /(/i)e-^V _1 d/i, A>0. (5.7) 
Jo 

Then for any d — 2s < (3 < d and any multiplication operator W > 

f°° f dr dt 

tT^F{W l l 2 B- p W 2 )< / eM-tBp){x,x)f(tW{x))—-. (5.8) 

Jo Jm d \ x \ t 

Note that the heat kernel exp(—tBp)(x,y) is well defined on the diagonal x = y by 
the semigroup property. Namely, exp(— tBp)(x, x) = J\exp(—tBp/2)(x,y)\ 2 \y\~^dy. 
For the proof of Proposition 15.31 one follows the proof of the CLR bound in |Llj (see 
also |Si2j and |RoSoj ) . As in the latter paper Trotter's product formula can be used 
in place of path integrals. For details we refer to Appendix |El 
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We shall now assume that F has the special form (|5.7jl in order to apply the trace 
estimate from Proposition 15.31 Given d — 2s < {3 < d and p = /?/(/? — d + 2s), 
Proposition 15.11 implies that 

ex${-tBp)(x,x)f{tWp(x)) 

< k m , s r [ r*f{tw fi {x)) lLr (U 



JR d \ x \^ t 

Note that Wp(x) p = K(x) p |x| /3 . We conclude that for any d/2s < p < oo, 

N(-l,H s -V)<K' pAs [ V(xf + dx, (5.9) 

where the constant is given by 

K' pAs = K M , S inf F(l)- 1 (j^r^mdt) . 

Here /3 = p(d — 2s) /(p — 1), and the infimum runs over all admissible functions / 
from Proposition 15.31 In order to obtain an explicit upper bound one may choose 
f(x) := (x — a) + and minimize over a > 0. 

Step 2. Now we use the idea of |LThj to deduce ()2.5|) from (|5.9|) . Fix 7 > and 
choose some d/(2s) < p < 7 + d/{2s). First we note that by scaling we have, for any 
t > 0, 

N(-r,H a -V)=N(-l,H a -V T ) 
where V T (x) := t~ 1 V(t~ 1 I 2s x). In view of ()5.9j) this yields 

N(-r,H s — V)< K' pds r-^ 2s [ V(x) p + dx. (5.10) 

Now, for any fixed < a < 1, one has by the variational principle 
N(-r, H s — V)< N{-{\ - a)r, H a — (V — ar)+) . 

Hence, by (jSHUJ), 

POO 

tr(H s - V)l = 7 / iV(-r, F s - ^)r 7 " 1 rfr 

< 7^,(1 - a)-^ 2s / / (y(z) - ar)l dxr^^' 1 dr . 

JO JR d 

We change the order of integration and calculate the r-integral first. For fixed x G M. d , 

/■oo 
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Here, B denotes the Beta-function B(a,b) = T(a)T(b)/T(a + b). Minimization over 
a £ (0, 1) and p £ (d/2s, 7 + d/2s) yields 



tr(H s - V)l < C d M [ V(x)l +d/2s dx (5.11) 

JR d 



with 



C djS ( 7 ):= , min L^K'^B^ + d/2s - p,p + 1) 

x (7 + d/2s - p)^" d/2s+p (p - d/2s)- p+d/2s \ . 
This concludes the proof of Theorem 12. II in the case A = 0. 



6. Extension to Magnetic Fields 

In this section we prove certain diamagnetic inequalities which allow us to extend 
the proof of Theorem 12.11 in the previous section to the case of non-zero magnetic 
fields. The main idea is contained in Proposition 16 .11 in Subsection 16. 11 The following 
Subsection 16.21 contains some technical refinements we will need. In Subsection 16.31 we 
describe the necessary modifications in the proof of Theorem 12.11 to include magnetic 
fields. The final Subsection 16.41 is devoted to the special case s — 1. 

Throughout this section we shall assume that A e L 1 2 oc (lR d ; IR d ) and that d > 2. Note 
that in d = 1 any magnetic vector potential can be removed by a gauge transformation. 

In Sj/3 = L 2 (M. d , \x\~Pdx) consider the quadratic form 

6/9^M:=/iM[kr ( ^ 2a)/2 «] + llkr ( ^ )/2 «l| 2 > veC™(R d \{0}). (6.1) 

In Appendix O we show that bp t A is closable and hence defines a self-adjoint operator 
Bp A in -ft/3- Our goal in this section is to show that exp(— tBp^) is an integral operator, 
whose kernel satisfies 

\exp(-tB/3 >A )(x,y)\ < exp(-tB p )(x, y) . (6.2) 

6.1. An initial inequality. We consider weighted magnetic operators u\D — A\ 2s uj 
where A e Lf oc (R d ) and u > with u + uj- 1 £ L°°(R d ). This is a self-adjoint operator 
in L 2 (IR d ) with form domain uo^ 1 dom \ D — A\ s . It satisfies the following diamagnetic 
inequality. 

Proposition 6.1 (Weighted diamagnetic inequality). Let d > 2 and < s < 1. 

Assume that A £ L 2 oc (R d ) and that u > with uo + uj' 1 £ L°°(lR d ). Then for all 
u £ L 2 (M. d ) and all t > one has 

I exp(-tuj\D - A\ 2s uo)u\ < exp(-tuj(-A) s uj)\u\. (6.3) 

Proof. First note that the assertion is true in the case u = 1, i.e., 

I exp(-t\D - A\ 2s )u\ < exp(-t(-A) s )\u\. (6.4) 
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Indeed, for s = 1 this inequality is proved in |Silj for all A e Lf oc (M. d ). The general case 
< s < 1 follows from the fact that the function A i— > e~ xs is completely monotone 
(i.e., its derivatives are alternating in sign) and hence is the Laplace transform of a 
positive measure by Bernstein's theorem |Doj . This reduces the problem to the case 
s = 1. 

Now assume that to is as in the proposition and write Ma := u\D — A\ 2s u. In view 
of the general relation 

exp(-tM A ) = s-lim (I + n~HM A )~ n (6.5) 

n— »oo 

it suffices to prove the inequality \(Ma + < (Mo + t) _1 |m|. But since (Ma + 

r)- 1 = u-\\D - A\ 2s + V)- 1 ^- 1 with V := tuj~ 2 it suffices to prove \(\D - A\ 2s + 
< (— A s + V) -1 |u|. In view of the relation 'inverse' to (|6.5|) . 

/>oo 

(|D- A| 2s + y)- x = / exp(-t(\D- A\ 2s + V))dt, 

Jo 

the assertion follows from (|6.4|) and Trotter's product formula. ■ 

Remark 6.2. We note that the diamagnetic inequality in the form ()6.4)1 implies (12.3)1 . 
This follows by integrating the square of ()6.4j) and evaluating the derivative with 
respect to t at t — 0. 

6.2. A refined inequality. In order to prove the desired inequality (|6.2j) . we have 
to extend Proposition 16. II in two directions. First, we want to use the singular weight 
u(x) = \x\ a , < a < s, and second, we want to replace the operator \D — A\ 2s by 
H SjA + /, i.e., we want to subtract the Hardy term. 

Recall that Bp t A was defined by ()6.1j) . The main result of this section is the following. 

Proposition 6.3 (Weighted diamagnetic inequality, second version). Letd > 2 
and d — 2s < f3 < d. Assume that A e L 2 oc (R d ). Then for all v G Sjp one has 

I exp(-t J B ( g >A )w| < exp(-iS j9 )|u|. (6.6) 

It follows from Proposition 15.11 that exp(— tBp) is an integral operator that maps 
L 1 (M d , Ixl'^dx) to L°°(M. d , Ixl^^dx). Hence ()6.6)1 implies that the same is true for 
exp(— tBp t A)- Moreover, the kernels are related by the inequality 1)6.2)1 . 

In the course of the proof we will need the following approximation result. 

Lemma 6.4. Let T n , T be closed, densely defined operators in a Hilbert space Sj with 
T n T* — > TT* in strong resolvent sense. Assume that there is a set T> C f"|domT n H 
dornT, dense in S), such that T n <p — > Tip for all <p £T>. Then T*T n —>■ T*T in strong 
resolvent sense. 

Proof. For 7 > and (p,ijj one has 

7 (<A (T n *T n + 7 )"V) = (<P, ^) ~ (T n ^ (T n T* n + 7 rXV0. 
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By assumption the right side converges to (ip,if)) — (Tip, (XT* + r y)~ 1 Tip) = j(<p, (T*T + 
7)~ 1 ?/>)> which proves that T*T n — > T*T in weak resolvent sense. However, the latter 
is the same as in strong resolvent sense, see jReSil[ Problem VIII. 20]. ■ 

Proof. Step 1. Consider again the unitary transformation U : L 2 (JsL d ) — > Sjp, which 
maps u I— > Ixl^ 2 ^. Then 

Qp, A :=U*Bp A U (6.7) 

is a self-adjoint operator in the unweighted space L 2 {R d ), whose quadratic form is 
given by 

qp )A [u] := h sA [\x\ a u] + \\\x\ a u\\ 2 , a = (/? + 2s - d)/2 . (6.8) 
Eq. ()6.6j) is equivalent to 

| exp(-tQ /3>A )M| < exp(-t<5 /3i o)|M|. (6.9) 

for u G L 2 (R d ). 

Step 2. We begin by considering the case where the 'potential terms' in the definition 
of Q/3,a are absent. More precisely, we consider the operator M^a in L 2 (R d ) generated 
by the quadratic form \\\D — A\ s \x\ a u\\ 2 on C^°(lR d \ {0}). We shall prove that for all 
u G L 2 (R d ) one has 

| exp(-tM PiA )u\ < exp(-tM A0 )|«|. (6.10) 
Let u n be a family of smooth positive functions which decrease monotonically to \x\ a 
and agree with this function outside a ball a radius n~ x . Similarly, for fixed n let uo n , m 
be a family of smooth positive and bounded functions which increase monotonically 
to uj n and agree with this function inside a ball a radius m. 

The operators \D — A\ s uo n and \D — A\ s uo ntm are easily seen to be closable on C^(R d \ 
{0}) (as in Appendix EJ) and we denote their closures by T n and T n>m , respectively. 
One finds that C °°(M d ) C domT n * with T*v = u n \D - A\ s v for v G C£°(R d ), and 
similarly for T* m . By construction of ui n ^ m the operators T n m T* m are monotonically 
increasing as m — > oo and hence converge in strong resolvent sense to T n T* by |ReSil| 
Thm. S.14]. Noting that T n>m <f — > T n (p for any cp G C^°{R d \ {0}) we conclude from 
Lemma E3I that T* m T n ^ m — > T*T n in strong resolvent sense. One checks that T^ m T n m 
coincides with the operator u) n>m \D — A\ 2s oj n ^ m from Subsection 16.11 and satisfies a 
diamagnetic inequality by Proposition 16.11 By the strong resolvent convergence the 
diamagnetic inequality is also valid for T*T n . Now we repeat the argument for n — > oo 
where we have monotone convergence from above. We apply [HI Thm. 1.2.3] noting 
that C^°(M. d \ {0}) is a form core for all the operators involved, and conclude that 
T*T n — > M a ^A in strong resolvent sense. This proves the diamagnetic inequality (|6.10|) . 

Step 3. Now we use another approximation argument to include the Hardy term. We 
define Rp,A via the quadratic form /i SiJ 4[|x| a u] on C^(R d \ {0}). Moreover, for n G N 
let W n (x) := C S: d mm{\x\~ 2 ( s ~ a \ n}. The boundedness of W n , ()6.10|) and Trotter's 
product formula show that the diamagnetic inequality is valid for Mg ^ — W n . Since 
Mp t A — W n — > R/3 t A in strong resolvent sense again by monotone convergence we find 
the diamagnetic inequality (|6.10j) with Mg ^ replaced by Rp^A- 
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Step 4- Finally, we note that Q@,a = Rp,A + \x\ 2a in the sense of quadratic forms. 
Moreover, C^°(M d \ {0}) is a core for both quadratic forms involved. Equation (jfi.9j) 
follows now from the diamagnetic inequality for -Rg,A by Kato's strong Trotter product 
formula [ReSil, Theorem S.21]. ■ 

6.3. Extension of Theorem 12.11 to magnetic fields. In the case of non-vanishing 
magnetic field, the proof of Theorem 12.11 is essentially identical to the one presented 
in the previous section. Although ()5.8|) does not necessarily hold with Bp^ instead of 
Bp, it does hold if Bp is replaced by -E>/3,a on the left side only! I.e., 

f°° C c\t dt 

tr^F(W^B^ A W^)< / / eM-tBp)(x,x)f(tW(x))^. 

JO JR d \ x \ 1 

For the proof, one uses the diamagnetic inequality ()6.2|) before applying Jensen's 
inequality (cf. Appendix |Ej). This leads to the conclusion that (j5.9J) also holds with 
magnetic fields, i.e., 



N{-1, H sA -V)< K' pAs [ V(x)l dx , (6.11) 

JM. d 



with the same (A-independent) constant K' p ds . 

For the remainder of the proof, we note that for any r > 0, 

N(-t, H sA -V) = N(-l, H sAt - V T ) 

where V T (x) := t^V^- 1 / 2 ^) and A T (x) := t- 1 ' 2s A(t~ 1 / 2s x). The scaling of A 
does not have any effect, however, since the constant in (j6.11j) is independent of A. 
Therefore (|5.11|) also holds with H s replaced by H s ^, with same constant Cd tS (l)- 

6.4. The special case s — 1. As noted in Remark 12.21 the proof of Theorem 12. II i ust 
given works also in the case s = 1 in dimensions d > 3. We briefly comment on the 
necessary modifications. 

The local Sobolev-Hardy inequalities for s = 1 have been proved in |B Vaj . Alterna- 



tively, one can obtain them following our proof in Section El Using the IMS formula 
()3.12|) one can obtain the global Sobolev-Hardy inequalities ()1.4|) . The rest of the 
proof goes through without change. To verify the Beurling-Deny criteria, one uses 
the ground-state representation ()4.1)1 instead of Proposition 14.11 Note also that the 
weighted diamagnetic inequalities in this section include the case s = 1. 

Appendix A. A Constant in the Sobolev Inequality ()2.6|) 

In this appendix we shall derive an explicit bound on the constants C q ^,\/2 for the 
Sobolev-Hardy inequalities ()2.6)1 in the case d = 3 and s — 1/2, which is of interest for 
our theorem on stability of matter. Let 3/2 < q < 3 and a := 2 — 3/q. For A > 1, let 

, . s 1 — a f°° dr . . 

P(Aj: "^A^A r (i+«)/2( r _A-2)- (A ' lj 
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We will show that 

<W < £ (1 - «) (S/4-) 4 / 3 mf - f ( " Q) . (A.2) 

6( 1 A>1 (l$i/2, 3 (a)| -P(A)) + 

We remark that in this special case 

2 , /7r(l-3/g)\ 



2 A 
|$i/2, 3 («)| = - - (1-3/g) cot ( - 



2 ) 

The estimate ()A.2|) is a consequence of the following two facts. First, we claim that 



,vr 2 



limsuplimsup/i 1/2 [V'A 5 ] < A 2(1 ~ q) — (1 - a) , (A.3) 

where ^ ,<5 (x) = A" -^'^/^) an d V ;£ ' <5 is defined for e, 5 > by 

{|x|~ a for \x\ < 1 , 

jxh 1 (1 -^(Ixl 2 - l) 5 ) for 1 < |x| 2 < 1 + l/e, (A.4) 
for |x| 2 > 1 + l/e. 

Note that ^ e,<5 does not satisfy the smoothness assumption of Lemma 13. 3[ but it can 
be approximated by such functions in /ii/2-norm. 

Secondly, we claim that p e ' s in (j3.10|) defined with the function if>\ — satisfies 

lim/' 5 (A) = p(A) (A.5) 

uniformly in 5 > and A > 1, with p(A) as in (jA.lJ) . Eq. ()A.2|) then follows from 
these two facts, proceeding as in the proof of Theorem 12 .31 Instead of choosing A such 
that p(A) = 1^*1/2,3(0!) |/2 we optimize now over the choice of A. 

For the proof of (|A.3j) we consider first an arbitrary radial function if> and, with a 
slight abuse of notation, we write if)(x) = ip{r) for r = Using the ground state 
representation in Proposition 14. 1[ introducing spherical coordinates, and integrating 
over the angles, we have 

r°° r°° rs 

hipty] = 8 / dr ds— — \rif)(r) - sif){s)\ 2 . 

Jo Jo [r — s ) 

By changing variables r 2 — > r and s 2 — > s, this yields 



h 1/2 [if)\ = 2 I dr I ds . IvW'(V^) " v^v^f • (A.6) 



Now assume that if) = ip e ' s as in (|A.4[) . By scaling, if suffices to prove (jA.3|) for 
A = 1. We split the integrals in (|A.6|) into several parts. First of all, we have 



dr I rfs ^^( r d-)/2_ s (i-«)/2) 2 = 2 f 1 ds L f 1 dt _}_( 1 _ t ^)/2f 
1 (r - s) 2 V ' J s<* J (1 - 1) 2 V ; 

(A.7) 
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This identity can be obtained by noting that the integral on the left is the same as 
twice the integral over the region r < s, and then writing r = st for < t < 1. Simple 
computations then lead to 

(i-t) 2{ > 4 J (i-tyj t J t 



(1 - a) /' ds— L^ 1 -*)/ 2 (1 - s ( 1 -")/ 2 ) . (A.8) 
Jo I — s 



We introduce the function 



We note that 77(A) = r"(A)/T(A) is the Digamma-function. It is then easy to see that 

1 ds-^—s {1 ~ a)/2 (1 - s^V 2 ) = v {2 -a)- 77(3/2 - a/2) . (A.9) 
1 1 s 



Altogether, we conclude that the contribution of r < 1 and s < 1 to the integral in 
u6]) is given by 



4 (77(2 -a) -77(3/2 -a/2)) . 
Similarly, we proceed with the other terms. We have 

lim / dr / cfe— — -^U 1 -^' 2 - \l -e 5 (s-l) 6 ]J 
£ Wo Ji (r-s) 2 V J +/ 

rfr / cfe^-i— (r^/ 2 - l) 2 = / dr— (r^ 2 - l) 2 . (A.10) 
- s) J Q 1 — r 



r — 



Here we have used dominated convergence, noting that the integrand is bounded from 
above by the L 1 function (r - s)~ 2 (2(l - r^ 1 "")/ 2 ) 2 + 2min{l, (s - 1) 2<5 }) for e < 1. 
The contribution of this term to (jA.6|) (noting that it appears twice) is thus given by 

4 (27/(3/2 -a/2) -77(1) -r/(2 -a)) . 

We are left with calculating 



l+l/e /"l+l/e 1 

dr / - (e\r - I) 5 - £ 5 (s - l) 5 ) 

1 7i (r - s) 2 



1 f 1 1 
rfr / ds (r 5 -s 5 ) 2 = 2 (77(1 + 25) - 7/(1 + 5)) . 



Jo 
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The last equality follows by proceeding as in (|A.7Jl - (|A.9Jl . The last term to evaluate 



is 



l+l/e rco ^ 



1 J l+l/e 



dr I ds X (l-e 5 (r-l) 5 ) 2 



r — s 



= / dr- (l - r s ) 2 = 2r](l + S) -7/(1) -77(1 + 25). 

JO 1 — T 

We have thus shown that 

lim/i 1/2 [^' 5 ] = 4 (^(3/2 - a/2) + 77(1 + 5) - 277(1)) . 

Using concavity of rj, together with t/(1) = 7r 2 /6, yields the estimate 



We proceed similarly for the calculation of p £ ' s . We have 

2 /"+'/' r"-)/' / l-^(r-l)^ , 

"^ = ^1 (T^lF^l 1 ^V^J*' 

Eq. (|A.5J) then follows by dominated convergence and integration by parts. 

Appendix B. Localization Formula in the Magnetic Case 

In this appendix we establish the analogue of Proposition 13.51 in the general case 
A 7^ 0. As explained in Section l2~2*l this is needed for the proof of Theorem 12.71 First 
recall that for s = 1 and Y^j=i Xj = 1 one nas 

/ \(D-A)u\ 2 dx = J2 \(D - A)( Xj u)\ 2 dx - ^TlVx/M 2 ^. 

JR d j=Q JR d JR d j=0 

In this case the localization error Y^j=o I^X?| 2 is local and independent of A. The 
analogue for s < 1 is 

Lemma B.l. Let d>2,0<s<l and A G L 2 oc (W d ). Then there exists a function 
kA on W 1 x R d such that the following holds. If Xo, ■ ■ ■ , Xn are Lipschitz continuous 
functions on R d satisfying YTj=oX 2 = 1; then one has 

n 

\\\D - A\ s u\\ 2 = ^2\\\D - A\ s Xj u\\ 2 - {u,L A u), u G dom \D — A| s , (B.l) 
i=o 

where La is the bounded operator with integral kernel 

n 

L A (x, y) := k A (x, y) ^(Xj( x ) ~ Xj(y)) 2 - 
3=0 

Moreover, for a.e. x,y G M. d 

\k A (x,y)\ < a s4 \x - y\~ d ~ 2s , and hence \L A (x,y)\ < L(x,y) 
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with L defined in Lemma \8. 51 

Proof. By the argument of |Sil| we can choose a form core for \D — A\ 2s which is 
invariant under multiplication by Lipschitz continuous functions. It suffices to prove 
(IB.lj) only for functions u from such a core. 
We write kA{x,y,t) := exp(— t\D — A\ 2s )(x, y) for the heat kernel and find 

n 

J2(XjU, (1 - exp(-t\D - A\ 2s )) Xj u) = (it, (1 - exp(-t\D - A\ 2s ))u) 



3=0 



\i>2jj k A(x,y,t)(xj(x) - X]{y)) 2 u{x)u{y) dxdy. 



+ 2 

3=0 

Now we divide by t and note that by our assumption on u the left side converges 
to ^j =0 |||-D — A| s Xj-u|| 2 as t — > 0. Similarly the first term on the right side di- 
vided by t converges to \\\D — A\ s u\\ 2 . Hence the last term divided by t converges to 
some limit (u, L^u). The diamagnetic inequality ()6.4j) yields the bound \kA(x,y,t)\ < 
exp(—t(—A) 2s )(x,y). This implies in particular that La is a bounded operator. Now 
it is easy to check that La is an integral operator and that the absolute value of its 
kernel is bounded pointwise by the one of L. ■ 

The following helps to clarify the role of the kernel fc^. 
Corollary B.2. Let u e dom \D — A\ 2s and assume that Q := IR d \ supp-u ^ 0. Then 

(\D - A\ 2s u) (x) = - / k A (x,y)u(y)dy forxeQ. 

Proof. Let cp G C£°(Q) and choose Xo,Xi suc h that xo = 1 on suppii, xi = 1 on 
supp<y9 andxo + Xi = 1- By polarization, (jB.lj) implies (<p, \D — A\ 2s u) = —(lp,Lau) = 
— J (p(x)kA{x, y)u{y) dxdy, whence the assertion. ■ 

Appendix C. Closability of the Quadratic Form bp^ 

This appendix contains some technical details concerning the quadratic form bp t A 
defined in (jfj.ljl . In particular, we shall show its closability. Throughout this appendix 
we assume that d>2, 0<s<l and A G L 2 oc (R d ). 

Lemma C.l. The sets C^(R d \ {0}) and T> := {w G dom(D - A) 2 n L°°(M d ) : 
suppw compact in IR d \ {0}} are form cores for \D — A\ 2s . 



Proof. It suffices to prove the statement for s — 1. In this case it is proved in jSilj 
that C£°(R d ) and V* := {w G dom(D - A) 2 n L°°(R d ) : suppw compact} are form 
cores for (D — A) 2 . Hence the statement will follow if we can approximate every 
function in any of these two spaces by functions from the same space vanishing in a 
neighborhood of the origin. But for functions u from C^°(lR d ) or T>* both functions 
Du and Au are square-integrable. This reduces the lemma to the case A = where it 
is well-known. ■ 
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Now let d — 2s < (3 < d and recall that the quadratic form qp t A was defined in 
flfi.Sjl . Note that C£°(R d \ {0}) is invariant under the unitary transformation in (j(17)) . 
Therefore, closability of qp^A and bp t A on C^°{R d \ {0}) are equivalent. 

Lemma C.2. The quadratic form qp^A, defined in \b\ <SJ) . is closable on C^°(M. d \ {0}). 

Proof. It suffices to show closability of the form r^jjit] := /i S) A[|a?| a w] on C^°(M. d \ {0}) 
for < a = (f3 + 2s-d)/2 < s. 

Let D be as in Lemma Kill We shall show that the quadratic form rp^A on T> is 
closable and that C^°(lR d \ {0}) is dense with respect to (r^M + IMI 2 ) 1 ^ 2 i n the 
closure of T> with respect to this norm. This implies the assertion. 

Let w G T>. Since \x\ a is smooth on supp«; and dom(D — A) 2 is invariant under 
multiplication by smooth functions we have \x\ a w G dom(D — A) 2 and hence \x\ a w G 
dom|.D — A\ 2s . Since \x\ a w has compact support it follows from Lemma lL3.ll and 
Corollary IB . 21 that |(|Z? — A\ 2s \ ■ \ a w){x)\ < C (w)\x\~ d ~ 2s for all large In particular, 
\x\ a \D-A\ 2s \x\ a w G L 2 (R d ). Moreover, \x\~ 2{s - a) w G L 2 (R d ). It follows that if u n G V 
such that u n — > in L 2 (M. d ), then the bilinear form associated with rp t A satisfies 

rp, A [un,v) = (un, \x\ a \D - A\ 2s \x\ a w - C s4 \x\- 2{s - a) w) -> 

as n — > oo. By standard arguments, this proves that is closable on P. 

In order to show the density of C^°(R d \ {0}) we let again w G V. Since \x\ a w G 
dom |Z?— A\ s , Lemma RTTI vields a sequence w n G C^°(lR d \{0}) such that || \D— A\ s (u n — 
\x\ a w)\\ + ||w n — — > 0. Hence, if w n := |x|~ a w n , then w n — > w in L 2 (lR d \ £>) 
and < r a ^A[w n — w] < \\\D — A\ s \x\ a (w n — w)\\ 2 — > 0. Moreover, Hardy's inequality 
implies also that || \x\~ s+a (w n — w)\\ — > and hence w n — > w in L 2 (B). This proves the 
density of C$°(M. d \ {0}) in the closure of V with respect to (V^aM + ||u>|| 2 ) 1/2 . ■ 



Appendix D. Nash's Argument 

For the sake of completeness we recall here how the Nash inequality (J5.4|) and 
the contraction property of exp(— tBp) imply the heat kernel estimate (|5.Hjl . Let 
k{t) := || exp(-tB p )v\\lp. Then 

k'(t) = -2b p [exp(-tB p )v] 

< -2C'j t M) l+1/p \\ zM-tBp)v\\-J^ dx) 

^ ~ 2 ^1 s M^) 1+1/p IMIli(^|-<W)> 
where we used (15 .4j) and the contraction property. Hence 

(Mt)- 1/p )'>2 P - 1 ^j|,|| L ^ | _^ ) 

and, after integration, 

Ktr 1/P > 2P- I ^ll^llx 1 t|-^)- 
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This means that 

\\ex V (-tB p )v\\l < (pC' q4 j2)n^v\\l 1(]xl .^ 
By duality, noting that Bp is self-adjoint in L 2 (\x\~@dx), this also implies 

|| expH^HI^,^ < (pC' qA j2)n-v\\v\\% p . 
Finally, by the semi-group property exp(— tBp) = exp(—tBp/2) exp(— tBp/2) 

|| expH^Hi^,-^ < (^W^^r^NI^,,-^). 
This is exactly the estimate (j5.3|) with the constant given in Proposition 15. 11 

Appendix E. A Trace Estimate 

In this appendix, we sketch the argument leading to Proposition ESI We emphasize 
that we shall ignore several technical details. We assume that W is smooth with 
compact support in R d \ {0}, and we put k(x, y, t) := exp(— tBp)(x, y). We claim that 
the trace formula 

tr % F(W 1/2 B^W 1/2 ) 

dt f f T~r , ( t \ / t , \ dxi dx n (E.l) 



Jo t n^ooJ Rd J^fJ- \ J J nj \nj^ 



3=1 v 7 \ 3=1 

holds true for any non-negative, lower semi-continuous function / vanishing near the 
origin. Here, F is related to / as in (|5.7|) . In the integral in (|E.1|) we use the convention 
that Xo = x n . Indeed, by an approximation argument it suffices to prove this formula 
for 

F(A)=A/(l + aA) , f(jx)=ne-°», 

where a > is a constant. Using the resolvent identity and Trotter's product formula, 
one easily verifies that 

F(W 1/2 B^W 1/2 ) = W 1/2 (Bp + aW)- l W l/2 

poo 

W l/2 exp(-t(B + aW))W l/2 dt 

lim T n (t) dt 
Jo n ~*°° 

in this case. Here, 

T n (t) := W 1/2 (exp(-tB /n) exp(-taW/n)) n W 1/2 . 

The latter is an integral operator. We evaluate its trace via integrating its kernel on 
the diagonal. Then we arrive at 



tr % T n (t) = J ■■■ J\{k^ 



nj \xi\p \x n \ 
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After symmetrization with respect to the variables this leads to 

1 f f tt , / t\ „ / t „ T , x \ dx-i dx 



tr % T n (t) = i j "... J f[k[ 




n 



Xj , Xj^i, 

n I \ n * — ' " I \x-\ r \x r . 

The claimed formula (|E.1|) follows if one interchanges the trace with the ^-integration 
and the n-limit. 

Now we assume in addition that / is convex. Then Jensen's inequality yields 

a j \ 11 — ' / \@ |x n |^ 

dx\ d/CCfi 

I r 1 \ P I t, I P 




j j=i v 



(Eq. ()E.1|) holds also in the magnetic case discussed in Section |Hl Before applying 
Jensen's inequality, one first has to use the diamagnetic inequality (|fi.2j) to eliminate 
the magnetic field in the kernel k, however.) Finally, we use the semigroup property 
to integrate with respect to the variables X2, ■ ■ ■ ,x n . We find that the latter integral 
is equal to 

dx 

k{x,X,t)f{tW{x))rrp. 

Plugging this into (|E.1J) leads to the estimate (j5.8j) . 

Details concerning the justification of the above manipulations can be found in 
RoSqj (see Theorem 2.4 there). 
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